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Abstract 

We show by an almost elementary calculation that the ADM mass 
of an asymptotically flat space can be computed as a limit involving 
a rate of change of area of a closed 2-surface. The result is essentially 
the same as that given by Brown and York [2] . We will prove this 
result in two ways, first by direct calculation from the original formula 
as given by Arnowitt, Deser and Misner and second as a corollary of 
an earlier result by Brewin for the case of simplicial spaces. 



1 Introduction 



It is well known that the for the Schwarzschild spacetime the ADM and the 
Schwarzschild masses are equal. Interestingly one can also easily recover the 
Schwarzschild mass by a simple calculation involving the areas of 2-spheres. 
Consider the time symmetric initial data for the Schwarzschild spacetime in 
the standard Schwarzschild coordinates 

ds 2 =(l- —) 1 dr 2 + r 2 dQ 2 
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The area of a sphere of radius r is A = 4nr 2 and its rate of change with 
respect to proper radial distance is 



dA 
ds 



( 2m 
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\ r 



) 



(1/2) 



For r » 2m we can expand this in powers of m/r to obtain 




The leading term on the right hand side is exactly the flat space value for 
dA/ds and thus we can re-write the above in the suggestive form 



with the subscripts 'f ' and 'c' denoting flat and curved space quantities re- 
spectively. The main result of this paper is that a similar result applies to 
all asymptotically flat spacetimes. 

The first point to address is how we might compute the curved and flat space 
versions of dA/ds for a general 3-metric. The procedure we shall follow is 
quite simple to state. We start with a given space (A4,g) and then chose a 
large closed 2-surface S in M.. For later convenience we shall assume that the 
surface S is everywhere convex (i.e. that its intrinsic Gauss curvature is ev- 
erywhere non-negative). This is not a serious constraint and should be easily 
satisfied provided the surface is built sufficiently far into the asymptotically 
flat regime of the space. 

The surface 5* will inherit, from its embedding in (A4,g), a unit normal n M 
and an intrinsic 2- metric h^. Thus we can compute not only the surface 
area A of 5 but also the surface areas for images of S obtained by dragging 
S a short distance along the integral curves of n M . This in turn will allow us 
to compute the rate of change of the area, dA/ds, with respect to the proper 
distance s off S. 

Now we shall engage in a small piece of surgery. First we discard the portion 
of (M.,g) outside of S (though we retain S itself). Since, by assumption, S 
is everywhere convex we are assured by the Weyl-Lewy embedding theorem 
[3J El that the truncated space may be extended by attaching a sub-set of 
flat space to S. The result is a new space which we denote by (Ai,g, S). 
The surface 5* can now be viewed as a transition surface from the curved 
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interior to a flat exterior. We shall use a '+' subscript to denote quantities 
on the exterior and likewise a '— ' subscript for quantities associated with the 
interior. Note that the surface components of the metric (i.e. the projection 
of g onto S) will be at least C° but no such continuity can be expected for 
the normal components (i.e. the projection onto the normal to S). Thus we 
can not expect the area A to vary smoothly as we cross S. Indeed we now 
claim that, in the limit S — > i°, the jump in dA/ds is directly proportional 
to the ADM mass, 

M =±i im (( d 4) -( d 4) 



8tt s^i° \ \ds ) , \ds ) 



We will present two very simple yet independent proofs. The first applies 
to smooth metrics and begins with the standard expression for the ADM 
mass in asymptotic Euclidian coordinates El Ej while the second proof, 
applicable to simplicial spaces, follows as a simple corollary to an earlier 
result by Brewin jH]. 

Brown and York PJ E] have produced a very similar result (with a slight 
change of notation), 

M = -L lim / (K. - K+) dS (2) 

where K + and K_ are the traces of the intrinsic curvatures of S+ and ir- 
respectively. By noting that for any 2-surface embedded in a larger ambient 
3-space, dA/ds = — J KdS it follows that two expressions for M are 
equivalent. 

What then is new in the present calculation? Simply that we offer an alter- 
native, and perhaps simpler proof of the main result (|T|). This fits the maxim 
that simple results demand simple proofs. 

Our calculations are adapted to just one quantity, the ADM mass. In con- 
trast Brown and York's computations are much more extensive. They derive 
quasi-local estimates for the energy and momentum and they show that these 
reduce to the familiar ADM expressions at spacelike infinity. Their calcu- 
lations are much more sophisticated than ours. In our approach we will 
compute the ADM mass directly from the standard integral formula. 
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2 The ADM mass 



Given a 3-manifold M. with an asymptotically flat positive definite 3-metric 
g the ADM mass of the space (Ad, g) can be defined by (0 El El) 



where S is a topological 2-sphere with outward pointing unit normal n M and 
area element dS. Note that the integrand is not in covariant from (indeed 
there is no covariant expression involving the first derivatives of the metric 
other than g^a = 0). Thus the above integral is defined only for a restricted 
class of coordinate systems and these are known as asymptotic Euclidian 
coordinates. In such a system the metric components are required to be of 
the form = 8^ + 0(1/ r). 

We first prepare for a series of local calculations by sub-dividing S into a finite 
set of non-overlapping patches Si, i = 1, 2, • ■ ■ . We do this simply to avoid 
any complications that might arise from doing a single global calculation. 
For much of the following we will focus on a single patch which we shall 
denote by Si. Our next step will be to show that a simple coordinate system 
can be constructed in a subset of A4 containing Si. The construction is quite 
elementary. Consider a typical Si of S and build a short cylinder Gj, with 
cross-section Si, by following the integral curves of the normal vectors to 
S. The cylinder will extend partly into A4.+ and partly into A4-. On the 
patch Si construct a 2-dimensional set of coordinates (u,v). If the cylinder 
is sufficiently short then each point P in the cylinder will be threaded by 
exactly one of the integral curves originating from Si. This allows us to 
assign coordinates (u, v, n) to P in the obvious way, namely n is the (signed) 
proper distance to P measured along the integral curve from Si and (u, v) 
are the coordinates in Si where the integral curve intersects Si. We take 
n > for points in Ai + and n < for points in We will now make 

a coordinate transformation in Ai + . Since g + is flat we know that we can 
transform the coordinates in Ai+ such that g + = diag(l, 1, 1). Finally we 
choose any smooth extension of the coordinate transformation into A4-. 
By this construction we will have built a (locally) asymptotically Euclidian 
coordinates in the cylinder Cj thus allowing us to use the above integral (jSJ) 
for the contribution from Si to the ADM mass M. 

Let Ai be the surface area of S^ Then using standard results from differential 
geometry we have 



M 



lim 

167T S^i° 




9^u, a ) n a dS 



(3) 
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where dS is the area element on Si and is the unit normal vector to Si. 
However 



= g^g^an" + <T (So.n^ - g^g^n" + g" u g, a n a !U 
= g^ (g^u, a - g av ,») n a + gTn^ + n% 
and thus we have 

J_dA = J_f \ <x dS+ ±_ , (g^n + n n ) dS 

8tt ds 16vr J s y vyM ' y 16vr ,/ s vy ' 



The integrand has been deliberately re-arranged to isolate the terms that 
arise in the original definition of the ADM mass. 

Let us now define [Q] for any object Q to be the jump in Q across S, that is 

[Q] = Q + -Q- 
Now consider the pair of integrals 

h = f [g^(g^«-g*^)n a ]dS 

J Si 

h = I [g^n^ + n^^dS 

J Si 

It is not hard to see that our claim (equation [TJ is true provided we can 
show that the second integral I 2 vanishes in the limit S — > i°. To prove 
this we need to carefully examine the asymptotic behaviour of the integrand. 
Let us first look at the relationship between and nt. Since n+ and nt 
are normal to the surface defined by n = we must have = f+n^ and 
n M _ = f-n^ where /+, /_ are scalar functions on Si. Thus we find 

where / = /_//+. We can also establish a similar result for the components 
n M . To this end, consider any pair of tangent vectors a M ,6 M to Si. Clearly 
= a^n+ and = b^n+ and thus n+ = q+e^dojbp where q + is a function on 
£». Likewise, for nt we find = q-e^i a a b/3 and thus = But 
since both normals are required to be unit vectors we must have q~/q+ = 1/f 
and thus 

n_ = ~J n + 
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What can we say about the asymptotic behaviour of /? We know that in the 
conformal analysis of asymptotically flat spacetime the 3-metric is C >0 at 
i°. The lack of differentiability is intimately tied to the direction dependent 
limits at i°. Thus it is not unreasonable to assume that the function / will 
vary smoothly along each spacelike geodesic ending at i°. The limit may well 
depend on which geodesic is chosen but along any one geodesic the function / 
should vary smoothly. Thus we propose that / should posses an asymptotic 
expansion of the form 



where a may be a direction dependent quantity near i . Alternatively one 
might argue that the following analysis is applicable only to those spaces for 
which / obeys the above form. In either case the class of spaces with this 
behaviour admits at least the time symmetric initial data of the Schwarzschild 
spacetime. 

Now let us turn our attention to the 3-metric and its behaviour across S. 
The integral curves through Si provide a natural means for foliating Cj by 
a sequence of images of Thus we can write the 3-metric in Cj in a 2+1 
form 



In this form h^ v is the induced 2-metric on each leaf of the foliation. Conti- 
nuity of the 3-metric across Si requires 



All indices will be raised and lowered using the 3-metric g^ u . 

We can now express all of the terms in I2 solely in terms of quantities based 
in Thus we obtain 




9fj,v h^ v -\- TljxTl, 



■v 



= [ VI 



i.e. h 






(1 



f)h> 
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where all terms on the right hand side are understood to be evaluated in Ai+ 
and thus, for simplicity, the + subscripts have been dropped. Note that we 
have also made use of the fact that the metric in A4+ is flat and thus in the 
Euclidian coordinates (which we chose earlier) h^ v n^^ = g^ v n^ v = n M iAt and 
= n^n u n^^. Combining the above we obtain 

brn^J + [<J = -J (1 - ff hTn^ 

Notice that the term is proportional to the trace of the intrinsic 

curvature on Si and thus varies as 0(l/r). Finally, since / — 1 = 0(l/r) we 
see that the integrand of 1% varies as 0(l/r 3 ) and thus I2 will vanish in the 
limit as Si — > i°. This completes the proof, that I2 vanishes. The last step is 
to form a sum over all of the Si,i = 1, 2, ■ • • yielding the formula as stated 
above (equation H}. 



3 Simplicial lattices 

Lattice theories such as the Regge Calculus (see [011113 E]) have enjoyed some 
success as a tool for numerical relativity though they have more commonly 
been used as a model for discrete theories of gravity. A popular example is 
the Regge calculus in which a simplicial lattice is built from a collection of 
non-intersecting blocks (which in 3-d are usually tetrahedral) , each with a 
flat internal metric, and an assignment of leg lengths to each leg in the lattice. 
Though the metric is flat in each tetrahedron, it is not flat everywhere. In 
fact the curvature must be treated as a distribution with support on the legs 
of the lattice. This complicates the study of such lattices but progress can 
be made. Indeed as an example of the use of standard distribution theory 
we showed [H] that the ADM mass of a simplicial lattice could be computed 



as 1 



M 



— V 9L (4) 



where the sum includes all the legs in the (triangulated) 2-surface S, 9 is the 
defect angle on the leg and L is the length of that leg. The defect angle 9 is 
defined as the 9 = 2n — £\ 9i where 9i are the angles subtended at the leg by 
the blocks attached to the leg. In our case we have four triangular cylinders 
attached to each leg on the surface S, with a pair of cylinders on either side 

1 Note that there was a typographical error in Brewin's paper |S] . The above corrected 
formula diffcr's from Brewin's by a factor of -1/2 
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of S. The defect angle can then be computed as 9 = —9 + — 9_ where 9 + and 
6L are the contributions from the respective pairs of cylinders in M. + and 
A4-. Finally, by inspection of Fig. 1, we can see that (9+)L = (dA/ds) + 
and (0-)L = (dA/ds)_. Upon substituting these into (jlj) we recover (JTJ). 

Equally we could have started with equation (JTJ) and used it to derive equation 

(ED- 
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This image contains two cylinders each with a triangular crossection. The 
outer cylinder, based on the vertices (a, b, c) and (a, b, c)+ is a cylinder of 
the simplicial lattice. The inner cylinder, based on the vertices (a, b, c) and 
(a, b, c)', is built from the simplicial lattice by dragging the lower triangle 
(a, b, c) outward along its normal a distance Ss. The areas of the two grey 
triangles are thus equal and thus the change in area from Si to Si+ equals 
the sum of the areas of the three (approximate) rectangles plus the three 
diamonds on S i+ . These later areas vary as 5s 2 and thus, to leading order, 
the contribution to 5 A from the leg (a, b) will be SA ab = 9 ab L ab 5s where 9 ab 
is the angle between the faces (a, b, a', b') and (a, b, a+, £>+). 
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